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PREFACE 

Since  the  establishment  at  the  turn  of  this  century 
chemical  engineering  has  steadily  kept  pace  with  the 
progress  in  other  branches  of  science  and  engineering. 
The  advances  in  molecular  physics  put  those  elegant 
molecular  concepts  at  the  disposal  of  chemical  engineers. 
Thus  they  are  able  to  relate  macroscopic  properties  of 
substances  to  microscopic  molecular  properties  to  obtain 
many  useful  correlations  and  better  design  procedures  that 
can  be  utilized  in  their  challenging  tasks.   The  statistical 
mechanics  enables  them  to  interpret  all  the  thermodynamic 
and  transport  properties  of  materials  in  terms  of  inter- 
molecular  forces.   Indeed  the  intermolecular  forces  play 
highly  important  roles  in  this  area.   Much  has  been  known 
but  much  has  yet  to  be  known  about  these  intermolecular  forces 

In  this  work  we  are  primarily  concerned  with  the 
prediction  of  properties  of  polar  and  nonpolar  gases  by 
means  of  intermolecular  potential  energy  functions.   In 
Chapters  I  and  II  the  hard-core  Morse  potential  is  investi- 
gated for  nonpolar  molecules  in  an  effort  to  find  a  potential 
model  which  can  be  applied  to  polyatomic  as  well  as  monatomic 
molecules.   In  Chapters  III  through  V  the  angle-averaged 
potential  model  with  the  temperature-dependent  parameters 

iii 


are  employed  to  correlate  properties  of  polar  gases  and 
mixtures  containing  polar  component.   The  parameters  of 
this  model  are  estimated  from  the  critical  temperatures 
and  polarizabilities .   In  Chapter  VII  transport  properties 
of  gases  are  predicted  from  saturated  liquid  densities 
through  the  principle  of  corresponding  states. 

The  reader  is  assumed  to  be  no  novice  on  the  subject. 
It  is  intended  that  each  chapter  is  complete  in  its  entirety, 
Therefore  no  exhaustive  review  is  included.   Instead,  each 
chapter  has  its  own  brief  introduction  at  the  beginning 
and  full  discussion  and  conclusion  in  the  course  of  the 
chapter.   However,  all  the  chapters  are  closely  related  to 
one  another. 

The  author  wishes  to  express  his  appreciation  to 
Dr.  K.  E.  Gubbins  for  serving  on  the  committee.  He  is 
especially  indebted  to  Dr.  T.  M.  Reed  who  directed  this 
research  and  served  as  the  chairman  of  the  supervisory 
committee.  He  is  also  responsible  for  influencing  the 
author  to  wander  in  the  jungle  of  "Molecular  Chemical 
Engineering  . 

Jae  Ho  Bae 
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This  work  is  concerned  with  the  pair-potential 
energy  functions  for  the  prediction  of  thermodynamic  and 
transport  properties  of  polar  and  nonpolar  gases. 

The  defect  in  the  Morse  potential  energy  function 
of  the  finite  value  at  zero  separation  is  avoided  by  assigning 
spherical  cores  to  molecules.   Tables  for  the  classical 
second  virial  coefficient  on  this  model  are  given.   The 
potential  parameters  for  ten  nonpolar  molecules  are  determined 
from  the  experimental  second  virial  coefficients.   The  experi- 
mental data  are  reproduced  very  well  and  the  agreement  is 
comparable  with  other  potentials.   This  model  is  also  extended 
to  12  binary  gas  mixtures  for  which  the  second  virial  coeffi- 
cients are  predicted  quite  adequately. 

For  polar  molecules  the  transport  properties  of 
dilute  gases  are  reproduced  accurately  by  the  12-6  pair- 
potential  energy  function  with  temperature-dependent  para- 
meter-, e  *  and  r*.   These  parameters  are  defined  from  the 
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classical  electrostatic  point  dipole-dipole  and  dipole-induced 
dipole  interactions  by  statistical  averaging  over  the  relative 
orientations  of  the  molecules,  and  adding  nonpolar  attraction 
and  repulsion  terms.   The  model  is  also  applied  to  binary 
mixtures  of  a  polar  and  a  nonpolar  gas  with  parameters 
determined  from  the  experimental  viscosities.   With  simple 
combining  rules  for  unlike  interaction  the  calculations  of 
viscosity,  diffusion  and  thermal  diffusion  factor  show  very 
good  agreements  between  experimental  and  calculated  values. 
The  usual  12-6  Lennard- Jones  integral  tables  are  used  with 
the  temperature-dependent  parameters. 

The  parameters  are  also  determined  for  this  model 
from  the  second  virial  coefficients  of  polar  gases.   The 
validity  of  these  parameters  is  further  confirmed  by 
computing  the  Joule-Thomson  coefficients  at  zero  pressure 
for  several  molecules. 

The  quantum  mechanical  expression  of  Mavroyannis 
and  Stephen  for  dispersion  attraction  confirms  the  empirical 
values  obtained  for  the  nonpolar  or  dispersion  attraction 
constant  for  small  polar  and  nonpolar  molecules.   From  this 
expression  and  the  principle  of  corresponding  states  a 
method  of  estimating  the  parameters  of  above  model  is  devised. 
The  estimated  parameters  give  reasonable  results  in  the 
calculation  of  both  thermodynamic  and  transport  properties 
of  polar  gases. 
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The  transport  properties  of  nonpolar  gases  are 
correlated  with  the  saturated,  liquid  density  parameters 
through  the  principle  of  corresponding  states.   Thus 
viscosity  and  diffusion  coefficients  are  predicted  within 
experimental  error  for  simple  molecules. 
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CHAPTER  I 


HARD-CORE  MORSE  POTENTIAL  FOR  THE 
SECOND  VIRIAL  COEFFICIENTS 


Most  of  the  potential  models  available  for  nonpolar 
molecules  predict  the  second  virial  coefficient  quite 
accurately  over  a  reasonable  temperature  range.   However, 
the  accuracy  becomes  worse  as  the  molecules  become  large 
and  nonspherical .   In  an  effort  to  find  a  potential  which 
can  be  used  for  large  nonspherical  as  well  as  for  the 
spherical  molecules,  the  hard-core  Morse  potential  is 
investigated  in  this  chapter. 

The  empirical  Morse  potential,  originally  proposed 
for  the  chemical  bond  energies  by  Morse,  x'  was  successfully 
employed  for  the  second  virial  coefficient  by  Konowalow, 
et_al.  Saxena  and  Gambhir^4'  extended  the  model  to 
the  binary  mixtures.  Saxena  and  Bahethi,  5'  Gambhir  and 
Saxena^6'  and  Saran^7'  used  the  model  for  the  correlation 
of  transport  properties  of  dilute  gases. 

In  the  hard-core  Morse  potential  the  convex  core 
model  of  Kihara^8'  is  combined  with  the  Morse  potential. 
Thus  the  potential  energy  of  interaction  is  assumed  to 
depend  on  the  shortest  distance  between  core  surfaces  and 
to  be  the  functional  form  of  the  Morse  potential.   The 


convex  core  can  have  any  shape  but  only  the  spherical  core 
is  considered  here.   The  potential  is: 


l=2i|  (r"rm) 


-2exp 
r  =  2a 


*bi)  ("m*] 


(r-r    )  >,r>  2a  (l) 


where  a  is  the  core  radius ,  o  is  the  separation  at  zero 
potential  energy,  e  is  the  minimum  potential  energy  at  r 
and  r  is  the  separation  between  the  core  centers.   The 
constant  c  has  the  same  meaning  as  in  the  Morse  potential. 
In  reduced  form  we  have: 

U*(r*)  =  U/e  =J  exp[-2c(r*-rm*)]-2exp  [-c(r*-rm*)J  ,  r*  >  a*    {2] 

r*  =   a* 


where 


o-2a 


(3) 


a*  = 


2a 


a- 2a 


It  can  be  easily  shown  that  the  reduced  separation  at  the 
potential  minimum  is: 


r  *  =  1  +  a*  +  (l/c)ln2 
m  v  /  / 


(4; 


Eq.  (k)    shows  that  r  *  of  this  model  is  the  sum  of  r  *  of 
1  m  m 

the  Morse  potential  and  the  reduced  core  radius.   In  fact 
Eq.  (2)  is  In  exactly  the  same  form  as  that  of  the  Morse 


potential  except  it  is  displaced  by  the  distance  a*. 

According  to  classical  statistical  mechanics  the 
second  virial  coefficient  for  angle-independent  potentials 
is  written  as: 


B(T) 


=  2  tt  N  j 


1-e     )r2dr 


(5) 


where  N  is  the  Avogadro  number,  8  is  1/kT,  and  k  is  the 
Boltzmann  constant.   By  the  substitution  of  Eq.  (2)  into 
Eq.  (5)  and  partial  integration  we  get  B(t)  in  reduced 
form  as  follows: 


B*(T*,c,a*)  =  a*3  +  3  /  (l-e-U*/T* 


r*  dr* 


where 


B*  =  B/ 
T*  =  kT/  e 


a* 


|  irN(  o  -2a)3 
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To  perform  the  numerical  integration  the  Gauss- 
Mehler  quadrature^9'  was  employed.   It  is  given  by: 

.+1 
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n   ,  _      -] 
J=l     J 


(8) 


where 


a  .  =  cos 
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2n 
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This  formula  is  accurate  and  simple  to  use/10'11'   In 


order  to  facilitate  its  use  we  introduce  a  new  variable 


x  = 


1+a* 
1+r* 


(10) 


Then  Eq.  (6)  becomes: 


+1 


B*(T*,c,a*)  =  a*3  +  3  J     '  g(x)(l-x2)+  l/2(l-x2),": 


dx 


(11: 


where 


g(x)  =(l-e-UVT*)(l+a*-x) 


2  1+a* 
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Since  r*  or  x  represents  the  distance,  g(x)  can  be  regarded 
as  an  even  function.   Finally,  by  Eqs .  (8)  and  (9)  we  obtain: 

B*(T*,c,a*)  =  a*3  +  2JL  "   g(a.)sin  I^^lL.  J   (14) 

n  J     l2n     i 


The  calculations  were  carried  out  on  the  computer 
and  the  values  of  B*(T*,c,a*)  are  given  in  Appendix. 
For  most  cases  n  =  SO  was  sufficient  to  give  one  part  per 
thousand  accuracy  but  for  the  extreme  cases,  e.g.,  small 
T*,  large  c  and  large  a*,  n  up  to  200  was  used.   The 
accuracy  is  believed  to  be  far  better  than  one  part  per 
thousand  for  most  of  the  values.   The  computer  program 
was  checked  with  a*  =  o  in  which  case  the  model  reduces  to 


the  Morse  potential.   The  result  showed  a  complete  agreement 
with  the  existing  tables^ 2>    for  all  the  values  of  T*  and  c. 

Potential  parameters  for  ten  nonpolar  molecules 
were  determined  from  empirical  second  virial  coefficients. 
At  first  the  parameters  were  approximated  by  a  graphical 
method  originally  used  by  Lennard- Jones . ^ 12'      They  were 
then  further  refined  on  the  computer  by  the  least-squares 
fit.  13;   it  was  assumed  that  the  errors  associated  with 
second  virial  coefficients  can  be  described  by  a  Poisson 
distribution.   Thus  a  weight,  the  inverse  of  the  observed 
B(T)  value,  was  given  to  each  point.   The  results  are  given 
in  Table  1  where  the  four  potential  parameters,  number  of 
data  points  used,  temperature  range  of  the  data,  root-mean 
square  deviation  and  references  to  the  data  are  given. 

The  validity  of  the  parameters  depends  on  the 
temperature  range  and  accuracy  of  the  data.   With  the  four 
parameters  of  the  present  model  the  data  in  a  narrow 
temperature  range  can  be  fitted  to  several  sets  of  parameters 
Among  the  molecules  considered  here  the  reported  data  on 
C3He,  n-C4Hi0,  and  SPS  were  not  consistent.   Thus  some 
sets  of  data  were  arbitrarily  discarded.   The  fit  was 
made  with  those  data  to  which  references  are  given  in 
Table  1.   Also  a  few  data  points  were  arbitrarily  rejected 
when  they  differed  from  the  calculated  values  by  more  than 
three  times  the  rms  deviation.   They  are  indicated  in 
Table  1. 


The  intermolecular  attraction  potential  derived 
from  dispersion  theory^ l4'  is  a  power  series  in  r-1. 
Dymond,  Rigby  and  Smith^15''  obtained  an  empirical  two 
parameter  potential  as  a  polynomial  in  r"1  for  the  inert 
gas  molecules,  and  quasi-spherical  polyatomic  molecules. 
An  inconsistency  of  the  Dymond- Rigby- Smith  model  is  that 
the  parameters  obtained  for  the  twelve  molecules  examined 
by  those  authors  do  not  fit  the  requirements  of  the 
corresponding  states  theorem  with  respect  to  reduced 
critical  constants  for  such  a  potential  function.   The 
Morse  potential  is  implicitly  a  polynomial  in  r_1  with 
three  parameters. 

A   difficulty  with  the  Morse  potential  is  that  at 
r  =  o  the  potential  has  a  finite  value.   Thus  at  high 
temperatures  molecules  on  collision  pass  through  one  another. 
This  defect  is  overcome  by  the  introduction  of  an  impenetrable 
spherical  hard  core,  the  size  of  which  is  a  fourth  parameter, 
to  give  the  hard  core  Morse  potential. 

The  parameters  given  in  Table  1  are  reasonably 
consistent  among  themselves.   As  the  molecular  size  increases 
the  value  of  c  increases,  which  corresponds  to  the  narrowing 
of  the  reduced  potential  well.   The  values  of  e /k  are 
always  larger  than  those  obtained  for  other  potentials .  le" 1Y' 
Compared  with  those  for  other  hydrocarbons  considered,  the 
value  of  e /k  for  CH4  seems  too  small,  but  the  rms  deviation 
is  much  less  than  that  for  other  molecules.   Also  the  core 
sizes  of  the  molecules  are  somewhat  large  compared  with 
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those  of  Kihara  potential.   For  example,  Sherwood  and 
Prausnitz  16   obtained  the  core  sizes  of  0.1841  A   for  Ar 

o 

and  0.7309  A  for  C02  while  in  this  calculation  we  obtained 
0.6738  A  and  1.0366  A  for  Ar  and  C02,  respectively. 
Considering  the  structure  of  the  molecules  the  former  values 
seem  too  small.   Since  the  present  model  has  a  finite 
potential  energy  at  r  =  2a,  it  must  have  larger  values  of 
a  and  a  than  those  of  the  Kihara  potential  to  have  the  same 
repulsion  contribution  to  the  second  virial  coefficient. 
The  core  sizes  are  inconsistent  among  Ar,  Kr,  SP6  and 
neoC5Hi2 . 

The  rms  deviations  between  calculated  and  experimental 
second  virial  coefficient  are  comparable  with  those  of  other 
potentials  for  simple  molecules . ( ls)   However,  for  the  hydro- 
carbons this  model  gives  better  agreement  than  the  unmodified 
Morse  potential . ( ly)   Nevertheless,  the  second  virial  coeffi- 
cients alone  are  not  sufficient  to  determine  the  unique  values 
of  the  four  parameters.   As  with  other  potential  functions, ^ 1Q) 
It  was  possible  to  find  several  sets  of  quite  satisfactory 
parameters . 

With  four  adjustable  parameters,  the  model  is  very 
flexible  and  from  the  results  of  this  calculation  it  can 
be  concluded  with  fairness  that  the  model  can  be  applied  to 
the  nonpolar  molecules  with  good  accuracy.   For  more  thorough 
comparison  between  potential  models  the  calculation  of 
transport  properties  for  the  present  model  is  highly 
desirable. 


CHAPTER  II 


SECOND  VIRIAL  COEFFICIENT  OF  BINARY  GAS  MIXTURES 
ON  THE  HARD-CORE  MORSE  POTENTIAL 


The  hard-core  Morse  potential  of  Chapter  I  is 
extended  to  mixtures  of  gases  since  there  exist  many 
experimental  data  on  the  second  virial  coefficient  of 
binary  gas  mixtures. 

In  order  to  approximate  the  potential  parameters 
for  the  unlike-pair  interaction  from  the  parameters  for 
the  like-pair  interactions  combining  rules  for  the  parameters 
should  be  used.   The  hard-core  Morse  potential  can  be 
written  in  the  form: 


where 


U(r)  =  A  e"2br  -B  e"br  (15) 


A  =  eexp  ("2Crm  \  (16) 

1  o -2a  J 

B  =  2  eexp  f^M  (17) 

I  a -2a  J 

(18) 


-2a 


and  all  the  symbols  have  the  same  meaning  as  before.   From 
the  theoretical  considerations  of  Zener^19'  the  parameter 


b  for  an  unlike  interaction  can  be  obtained  by: 

bX2  =  -  (bn  +  b22)  (19) 

where  subscripts  denote  interactions  between  molecular 
species.   Unfortunately,  no  combining  rules  for  the 
parameters  I  and  B  can  be  obtained  from  theory.   Therefore, 
the  empirical  rules  of  Mason  and  Rice^  and  Saxena  and 
Gambhir  4   are  adopted  for  these  parameters.   Thus, 

-  -   1/2 

A12  =  (AUA22)  (20) 

-      -  _   x/2 

B12  =  (B11B22)  (21) 

Since  the  cores  are  impenetrable,  rigid,  hard  spheres, 
it  is  reasonable  to  use  the  arithmetic  mean  for  the 
parameter  a. 

ai2  =  -  (an  +  a22)  (22) 

2  v  ' 

Prom  Eqs.  (20)  and  (21)  it  follows: 

i/a 
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=  (  £u  £22)  (23) 


The  derivation  of  Eq.  (23)  is  rather  empirical  but 
Eq.  (23)  is  the  common  combining  rule  for  12-6  potential 
Following  Reed,(21>22)  Hudson  and  McCoubrey, ( 2s)  and 
Mavroyannis  and  Stephen^ 24^  we  have: 


1/ 


=  (  en  E22)"  "  fT  f6  (24) 

-L   s 
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where 


a   a  _  \   1/4 
f  x   =  2  ' — L- 

ZiZ2 


'[(=• 


(25) 


f   =  2(V 
s    s  c 


1/6  ,  r 
/Vc2)    /[ 


/e  .  r  1/3' 

1  +  (v  /v 


ci   c2' 


(26) 


and  a  is  polarizability,  V  is  critical  volume  and  subscripts 
denote  molecular  species.   From  Eq.  (19)  the  following 
relationship  is  obtained: 


C: 


i2-2ai2 


n-2a: 


22~2a22 


(27) 


Combining  Eqs .  (20)  and  (27)  with  the  expression  for 
r*   of  Eq.  (4)  we  obtain  the  combining  rule  for  a 12 
as  follows: 


Cn  o 


C22  g 22 


<m-2an       °  22-2a22  I V  °n-2a 


22-2a22 


(28; 


Consequently  Eqs.  (22),  (23)  or  (24),  (27)  and  (28)  give 
the  desired  combining  rules  for  the  parameters  &X2,      Ei2, 
C12  and  °i2. 

The  cross  second  virial  coefficients  are  calculated 
for  12  binary  gas  mixtures.   For  the  purpose  of  comparison 
calculations  were  made  on  the  Morse  potential  and  Lennard- 
Jones  12-6  potential  whenever  parameters  were  available. 
For  the  hard-core  Morse  potential  the  parameters  in 
Chapter  I  were  used  and  for  other  potentials  the  parameters 
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used  are  listed  in  Table  2.   For  the  hard-core  model  the 
suggested  combining  rules  are  used  and  for  the  unmodified 
Morse  potential  the  same  rules  are  used  with  a  =  o  in  which 
case  Eq.  (27)  reduces  to  arithmetic  mean  for  C12.    Saxena 
and  Gambhir's  combining  rule  for  Morse  potential  implies 
implicitly  the  arithmetic  mean  for  C12  even  though  they 
apparently  did  not  notice  it.   For  12-6  potential  model 
the  usual  geometric  mean  for  energy  parameter  and  arithmetic 
mean  for  size  parameters  are  used.   All  the  calculations 
were  made  with  and  without  f-f actors  for  e  12,  I.e., 
by  Eqs.  (23)  and  (2k).      The  critical  volumes  and  polariza- 
bilities  needed  for  f-f actors  are  given  in  Table  3.   The 
results  which  give  better  agreement  to  the  experiments 
are  summarized  in  Table  4.   The  cases  when  f-factors  are 
used  are  indicated  by  a  small  f  on  the  right  side  of  the 
first  row  of  that  system. 

For  all  12  systems  considered  here  the  overall  rms 
deviation  in  cc./mole  is  8.56  and  11.44  for  the  hard-core 
model  and  the  12-6  model,  respectively.   For  the  five 
systems  for  which  unmodified  Morse  potential  is  used,  the 
deviation  is  8.29,  5.13  and  9. 80  for  hard-core  model,  Morse 
and  12-6  potential,  respectively.   It  is  noted  from  this 
calculation  that  the  two  parameter  Lennard-Jones  potential 
is  inferior  to  three  parameter  Morse  potential  or  four 
parameter  hard-core  Morse  potential. 
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The  results  show  an  apparent  superiority  of  the 
unmodified  Morse  over  hard-core  Morse  potential  in  the 
mixture  calculations.   However,  with  a  close  look,  a 
definite  conclusion  cannot  be  drawn.   It  is  well  known 
that  the  second  virial  coefficient  is  not  sufficient  to 
determine  unique  values  of  several  parameters  of  a  potential 
model.   The  validity  of  the  parameters  depends  on  the 
temperature  range  and  the  accuracy  of  the  data.   It  should 
he  pointed  out  that  for  the  Morse  potential  all  the  tempera- 
tures in  the  mixture  calculations  lie  within  the  temperature 
range  for  which  parameters  were  determined  while  for  hard- 
core model  this  is  not  true. 

Among  the  systems  considered  here  Ar  +  Kr  is  the  only 
system  for  which  the  parameters  for  pure  components  were 
determined  from  the  same  data  for  both  models.   These  systems 
show  essentially  the  same  results.   For  CH4  +  Ar,  mixture 
calculations  were  made  from  108.6°K  to  295°K  while  CH4 
parameters  of  hard-core  Morse  potential  were  determined 
from  the  data  of  273°K  to  623°K.   At  295°K  the  calculated 
cross  second  virial  coefficient  on  the  hard-core  agreed 
completely  with  the  experiment.   This  trend  is  also  true 
with  CH4  +  C4H10  system.   In  the  temperature  range  of 
parameter  determination,  the  hard-core  model  gives  better 
results  than  the  Morse  potential  in  the  above  cases.   As 
for  C3H3 ,  it  showed  the  largest  deviation  in  the  parameter 
determination  for  the  hard-core  model  and  from  the  trend 
of  the  systematic  deviations  in  CH4  +  C3H8  system,  the 
validity  of  the  parameters  for  C3Ha  is  rather  suspected. 
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It  seems  that  the  combining  rules  are  satisfactory, 
but  it  is  interesting  that  for  one  potential  model,  when 
better  agreement  is  obtained  by  Eq.  (23)  than  by  Eq.  (2*0, 
the  same  is  true  for  other  potentials.   Even  the  empirical 
correlation  of  Huff  and  Reed^2-   followed  the  same  pattern, 
It  appears  that  regardless  of  the  functional  form  of  the 
potential  energy  function,  the  energy  of  potential  minimum 
for  unlike  interaction  can  be  approximated  from  those  for 
like  interactions  by  the  same  combining  rule.   However, 
neither  Eq.(23)  nor  Eq.  (24)  seems  to  present  the  unlike 
interaction  properly.   Actually,  the  dispersion  attraction 
contains  terms  of  higher  order  of  inverse  separation  than 
6.   Eq.  (23)  which  was  derived  on  r~6  term  only  is  not 
truly  theoretical.   To  account  for  the  unlike  interaction 
properly,  a  more  elaborate  combining  rule  for  the  energy 
parameter  should  be  developed. 

It  is  felt  that  for  more  critical  comparison  of 
potentials,  transport  properties  on  the  hard-core  Morse 
potential  should  be  calculated.   Until  that  time  the  final 
comparison  should  be  reserved. 


CHAPTER  III 


ANGLE- AVERAGED  PAIR- POTENTIAL  ENERGY  FUNCTION  FOR 
DILUTE  POLAR  GAS  PROPERTIES 


While  there  have  been  many  pair-potential  energy- 
models  for  spherical  nonpolar  molecules,  the  angle-dependent 
12-6-3  Stockmayer  potential  is  the  only  one  that  has  been 
used  extensively  for  the  interaction  of  a  pair  of  polar 
molecules.   Rowlinson^ ze'    demonstrated  the  successful 
correlation  of  second  virial  coefficients  of  polar  gases 
with  this  potential.   Saxena  and  Joshi^27^  extended  the 
model  by  using  18-6-3  and  28-7-3  powers  on  the  reciprocal 
of  the  intermolecular  separation.   O'Connell  and  Prausnitz^28 
modified  the  Stockmayer  potential  by  introducing  the  Kihara 
hard-core  model  for  the  nonpolar  part  of  the  interaction 
in  place  of  the  Lennard- Jones  model.   Monchick  and  Mason 
have  applied  the  12-6-3  Stockmayer  potential  to  the  correla- 
tion of  transport  properties  of  dilute  gases. 

A  potential  energy  function,  which  may  be  called  the 
12-6-6  model,  for  polar-polar  and  for  polar-nonpolar  pairs 
contains  the  statistically  weighted  mean  potential,  averaged 

over  relative  orientation,  for  the  classical  electrostatic 

( «i  ) 

dipole-dipole  and  dipole-induced  dipole  interaction^  x/ 

together  with  the  Lennard- Jones  12-6  potential  for  the 
contribution  of  the  dispersion  force  interaction: 


ill 


29, 3o) 
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<  <j)  >    =   e° 


U  [(r°i/^12  -*<*V*>S]  -[2*1**    3/3 

u  2  a  1/r6  (29) 

.1   iJ 


2      , 
V    .   a    + 

1   J 


<4»>1.  is  the  total  weighted  average  mutual  potential 
energy  of  the  pair  i,  j.  e  °   and  r° .  .  are  empirical 
constants  characteristic  of  the  dispersion  interaction 
of  the  molecular  pair  i  and  j,  p  is  dipole  moment  in  zero 
external  field,  a  is  mean  polarizability,  and  Bis  (kT)"1. 
It  is  assumed  in  the  derivation  of  the  polar  interaction 
term  that  molecules  are  point  dipoles  in  the  homogeneous 
field  produced  by  the  molecule  at  the  other.   The  effective 
potential  energy  of  the  pair  of  molecules  is  approximated 
by  the  term  in  r  s  in  the  orientation-averaged  electrostatic 
interaction. 

<<j>>  is  thus  an  angle-independent  potential  but  it 
is  a  function  of  temperature  as  required  by  the  statistical 
averaging  over  the  relative  orientations. 

The  right-hand  side  of  Eq.  (29)  may  be  arranged^31) 
into  the  usual  12-6  bireciprocal  form, 


<♦>  .  .  =   e*.  . 

IJ       IJ 


r*.  ,  lis 

-^M    -2  U 


(30 
where   e*±j  and  r*   are  functions  of  temperature, 
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and 


ij      ij 


14-  u  2y  2  6/3r°?.e  ° 


1   J 


1J    ij 


2,2 

+  (U    a+y    a)/2e°   ros 

i   J      J   i       ij    Ij 


r*.  .  =  r°.  .  (  e  °.  ./E  *.  .) 


i/i; 


(51: 


(32; 


This  potential  obviates  the  calculation  of  any  additional 
tables  of  reduced  virial  coefficients  and  collision  integrals 
for  dilute  gas  properties  since  such  tables  already  exist^32' 
for  the  functional  form  of  Eq.  (30).   In  this  model  the 
potential  is  averaged  over  ail  orientations  before  the 
collision  integral  is  calculated,  whereas  in  the  model  of 
Monchick  and  Mason^29^  the  collision  integral  is  calculated 
as  a  function  of  orientations  and  then  averaged  over  all 
orientations,  giving  equal  weight  to  each  orientation.   In 
both  models  it  is  assumed  that  there  is  no  interchange  of 
rotational  and  translational  energies  during  a  collision. 
The  Stockmayer  potential  carries  no  terms  for  the  dipole- 
induced  dipolar  interaction  that  is  included  in  the  potential 
of  Eq.  (29).  <*>is  a  simple  function  and  as  shown  here,  is 
as  good  as  or  better  than  the  Stockmayer  potential  for 
correlating,  interpolating,  and  extrapolating  dilute  gas 
data.   It  should  be  clearly  understood  here  that  the  present 
model  is  distinctly  different  from  the  12-6  model  in  that 
the  averaged  dipole-dipole  and  dipole-induced  dipole 
interactions  are  superposed  to  the  latter  to  obtain  the  former 
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For  a  single-component  system  1  =  j  always  and  the 
species  subscripts  may  be  deleted. 

To  calculate  properties  of  a  single-component  gas 
the  parameters  e  °  and  r°  as  well  as  the  dipole  moment 
and  the  (mean)  polarizability  are  required.   The  potential 
parameters  are  obtained  by  fitting  the  theoretical  relation- 
ships to  the  experimental  data.   It  is  assumed  that  the  dipole 
moment  and  polarizability  are  known  molecular  constants. 
The  data  used  in  the  present  calculations  are  those  of  gas 
viscosities  treated  by  Monchick  and  Mason.  29' 

The  rigorous  kinetic  theory^  J/  gives  the  following 
expession  for  first  order  approximation  to  viscosity, 
n in  gm./cm.  sec.,  in  a  dilute  gas: 

1/2 
nxl07  =  266.92  ^^  (33) 

a2  fi*(2'2) 
where  M  is  molecular  weight,  T  is  temperature  in  °K  , 

o 

0  is  the  separation  in  A  at  zero  potential,  and  fi*  is 
the  dimensionless  collision  integral.   Here  °  can  be 
replaced  by  r*,  the  separation  at  the  potential  minimum, 
using  the  relationship  2  -g  a    =   r*  for  the  12-6  potential. 
Since  the  experimental  error  is  In  the  same  order  of 
magnitude  as  that  of  the  correction  term  for  higher 
theoretical  approximations,  the  first  approximation  is 
adequate  here.   The  two  Eqs.  (31)  and  (32)  contain  four 
unknowns,  e*,  r*,  e°,  and  r° .   Choosing  two  temperatures, 
Ti  and  T2,  we  obtain  ri*  and  r2*  in  terms  of   ex*  and   e2* 
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from  these  relationships  as: 


ri' 


£l*   £2*)2 


(  *  x   -     Ba) 


(34: 


where  the  subscripts  1  and  2  now  refer  to  the  temperatures. 
The  same  relationship  for  r2*6  can  be  obtained  by  interchange 
of  subscripts  1  and  2.   The  r*'s  are  eliminated  from  these 
equations  by  introducing  the  theoretical  expression  from 
Eq.  (33)  for  each  temperature.   The  resulting  expressions 
rearranged  give  the  working  equations: 


*x*     e2*  =  (  ^*  -  Ai) 


15) 


E  ,  •*  E  „* 


E  2*  -  A2) 


^l2 


where 


7631x10  ^     P4  n*3  n  a 


n  i"  n  i 


(MTi)3//j 


1  -   p  2 


(36) 


)    (37) 


and  similarly  for  A2 . 

When  viscosities  at  Ti  and  T2  are  known  then  e  x* 
and  £  2*  are  the  only  explicit  unknowns  in  Eqs.  (35)  and 
(36)  .   The  values  of  %.*  and  e  2*,  which  satisfy  these 
two  are  the  energy  parameters  at  Ti  and  T2 .   ri*  and  r2* 
are  then  evaluated  by  Eq.  (3^) . 

The  temperature- independent  parameters,  r°  and  e ° , 
are  obtained  in  terms  of  r*  and  e*  from  Eqs.  (31)  and  (32)  as: 


1  - 


r*3 


*  p*6 


(38) 
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e   °  =   e*(r*/r0)12  (39) 

In  this  way   e°  can  be  computed  from  two  experimental 
viscosities . 

The  potential  parameters  for  several  polar  gases 
determined  by  this  method  are  given  in  Table  5,  along  with 
molecular  constants  y  and  a  ,  and  references  for  viscosity 

data.   Dipole  moments  are  taken  from  Monchick  and  Mason^2 

( 
and  polarizabilities  are  taken  from  Landolt-Bornstein  Tables;3 

Eqs.  (35)  and  (36)  each  gave  curves  for   ei*  vs.   e2* 

of  small  curvatures  and  a  small  portion  of  each  curve  was 

approximated  by  a  straight  line.   A  graphical  method  was 

used  to  find   ex*  and   £2*  In  this  work.   In  all  cases, 

except  HBr  for  which  there  are  only  two  viscosities,  the 

intersection  of  the  two  curves  could  be  easily  and  accurately 

located.   For  water  vapor   e*  is  very  large  at  low  temperatures 

(  e*/k  =  228o°K  at  273°K) .   To  be  able  to  use  the  12-6  tables 

for  water  vapor,  the  high  temperature  data  of  Bonilla,  Wang 

(3s) 

and  Weiner's  smoothed  datav  '    were  used  without  any  correc- 
tion mentioned  by  others. ^Z3> 

The  experimental  and  calculated  gas  viscosities  at 
one  atmosphere  pressure  are  compared  in  Figures  1,2  and  3- 
The  data  points  used  in  the  determination  of  the  parameters 
are  marked  by  crosses.   Arbitrarily  chosen  data  points  are 
reproduced  in  the  figures.   The  curves  are  calculated  values. 
For  most  of  the  molecules  considered  the  agreement  is  within 
the  experimental  error.   The  deviation  in  H2O  is  the  largest 
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among  all  the  molecules  studied.   At  the  extreme  high 
temperature  of  2073°K  the  deviation  is  +3.8^.   For  the 
remaining  molecules  the  agreement  seems  to  be  as  good  as 
those  obtained  for  nonpolar  gases  with  the  usual  potentials /32^ 

The  method  of  determining  parameters  from  experimental 
data  is  straight-forward.   It  does  not  suffer  from  the 
subjectivity  frequently  characteristic  of  graphical  methods. 
Theoretically  E  *  and  r*  should  be  constant  at  a  given 
temperature.   We  paired  one  temperature  with  several  others 
and  calculated  e*   and  r*  at  the  fixed  temperature  from  the 
several  different  combinations.   The  values  found  differed 
less  than  1%   from  one  another. 

The  temperature-dependent  factors  in  Eos.  (31)  and 
(32)  are  monotonously  decreasing  and  increasing  respectively. 
The  temperature  dependence  is  rather  small  for  r*  and  large 
for  e*,    especially  when  the  molecule  has  large  dipole  moment. 
For  the  feebly  polar  molecules,  like  HI  and  CHC13,  E  */  e° 
is  near  unity  while  for  strongly  polar  molecules,  like 
H20,  NH3,  and  (CH3)2CO,   e*/  e°  is  large;  e.g.,  for  H20 
e  */  e°  =  11.4  at  273°K. 

The  second  virial  coefficient,  computed  with<<j>>   and 
the  12-6  tables  and  the  parameters  of  Table  5,  are  shown 
in  Table  6.   The  second  virial  coefficient  data  for  S02  are 
taken  from  Kang  et  al.(36>   and  all  others  from  Rowlinson. (2e) 
The  agreement  between  observed  and  calculated  values  are  as 
good  or  as  poor  as  those  obtained  by  Monchick  and  Mason^29^ 
with  the  angle-dependent  potential  and  viscosity  parameters. 
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No  potential  has  been  found  which  is  successful  with  a  single 
set  of  parameters  for  a  given  molecular  pair  for  both 
equilibrium  and  transport  properties  of  polar  molecules. 

Usually  the  transport  properties  are  better  predicted 
with  the  parameters  determined  from  other  transport  properties 
For  the  polar  molecules  under  consideration  the  only  available 
experimental  value  is  the  self-diffusion  coefficient  of  HC1 . 
Theoretically  the  self-diffusion  coefficient,  D  in  cm.2/sec, 
in  the  first  approximation  is  written  as: 


=  0.002628  (T3/M? 


i/a 


D  =  0.002628  \x    '  lV , (k0 


P  a2  n 


*U,i) 


where  P  is  the  pressure  in  atmosphere  and  other  symbols 
have  the  same  meaning  as  before.   By  use  of  this  formula, 
the  self-diffusion  coefficient  of  HC1  at  295°K  was  computed 
to  be  0.1261  from  12-6  tables  and  the  parameters  of  Table  5. 
If  we  take  Braune  and  Zehle's^37'  value  of  0.1254  (as 
corrected  for  the  mass  of  isotopic  tracers  by  Monchick 
and  Mason)  the  deviation  is  +  0.56$  which  is  well -within 
the  experimental  error. 

The  parameters,  z°   and  r°,  hypothetically  are  those 
for  the  dispersion  part  of  the  total  intermolecular  interactioi 
The  theory  of  Mavroyannis  and  Stephen^24'  for  interaction 
between  atoms  gives  an  expression  for  the  interaction  constant 

£°r°s  In  terms  of  the  polarizability,  a,  and  the  total 
number  of  electrons,  Z,  in  the  molecule  as: 


e   °r°S  =  6.292xl0"12  a  */z   Zl/s   erK  X  1' 


(41) 


when  12  is  used  as  the  repulsion- term  exponent.   The 
negative  of  four  times  E  °r°6  is  the  coefficient  on  r"6 
in  the  dispersion  energy  interaction.   This  theoretical 
value  for  e°r°°  is  compared  in  Table  7  with  the  corresponding 
empirical  value  calculated  with  the  parameters  of  Table  5. 
The  agreement  is  rather  close  for  the  smaller  polar  molecules 
and  becomes  poor  for  the  larger  polyatomic  ones. 

The  agreement  between  the  Mavr-     '•;- Stephen  theory 
and  data  for  nonpolar  molecules  is  illustrated  in  the  lower 
part  of  Table  7.   The  agreement  is  very  good  for  the  monatomic 
and  diatomic  molecules  and  becomes  poor  as  the  number  of 
atoms  per  molecules  increases. 

The  parameters  e °  and  r°  for  the  polar  molecules 
thus  appear  to  fulfill  their  role  as  representative  of  the 
dispersion  part  of  the  total  interaction  in  a  manner  which 
is  as  realistic  as  the  parameters  for  the  toal  interaction 
between  two  nonpolar  molecules  in  the  12-6  form.   It  is 
interesting  to  find  that  the  Mavroyannis-Stephen  expression 
for  dispersion  interaction  so  accurately  confirms  the 
experimental  interaction  constant  £°r°6  for  small  molecules, 
both  polar  and  nonpolar.   The  London  formulas,  on  the  other 
hand,  predict  values  of  £°r0°  that  are  low  by  a  factor  of 
1.6  for  monatomic  molecules  and  still  worse  for  diatomic 
and  larger  molecules.'24' 


CHAPTER  IV 


TRANSPORT  PROPERTIES  OF  POLAR- NONPOLAR 
GAS  MIXTURES 


The  angle-averaged  potential  energy  function  of 
Chapter  III  Is  applied  to  binary  mixtures  containing  polar 
gases,  specifically  to  polar-nonpolar  mixtures  since  no 
data  exist  for  polar-polar  mixtures. 

The  unlike-pair  interactions  may  be  approximated  from 

those  for  like-pairs.   It  is  clear  from  Eq.  (31)  that  e  *. . 

cannot  be  a  geometric  mean  of  e  *..  and  e  *...   The  dispersion 

11        JJ 

force  parameter  e  °. .  is  expected  to  be  close  to  geometric 

mean  of  e ° . .  and  e0..  and  the  arithmetic  mean  rule  is  used 
11        JJ 

for  r° 

ij 

r°.  .  -  (r°..  +  r°  .  .) /2  (45) 

Using  dipole  moments  and  polarizabilities,  e  *in-  and  r*±i 

can  be  calculated  by  Eqs .  (31)  and  (32). 

When  one  member  of  the  pair  is  nonpolar,  e  *. .  and 

r*   are  no  longer  temperature-dependent  since  the  t 
ij 

p  2  y  2  is  then  zero.   Nevertheless,  e  *.  .  is  not  a  geometric 

i    J  XJ 

mean  since  the  dipole-induced  dipole  interaction  is  included 
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in  the  potential.   Usually  the  induced  moment  is  ten  times 
smaller  than  the  inducing  one  and  the  model  is  almost  the 
same  as  the  12-6  potential.   However,  this  should  be  more 
realistic  than  the  12-6  model  and  as  shown  here  the 
parameters  for  the  12-6-6  model  successfully  predict  the 
properties  of  polar-nonpolar  mixtures. 

Making  use  of  the  combining  rules  in  Eqs.  (42)  and 
(43)  and  the  empirical  values,  previously  obtained  from 
gas  viscosities,  of  E°  and  r°  in  like-pair  interactions, 
the  viscosity,  diffusion  coefficient,  and  thermal  diffusion 
factor  were  calculated  and  compared  with  experimental 
properties  for  the  binary  systems  in  Tables  9,  10,  and 
Figures  4,5,  and  6.   The  potential  parameters  for  nonpolar 
gases  used  in  the  calculation  are  given  in  Table  8  along 
with  polarizabilities.   Thermal  conductivity  is  not  considered 
here  because  of  its  dependence  on  the  internal  energy  transfer, 
From  the  kinetic  theory  the  expression  for  the  viscosity 
of  binary  mixture  can  be  written  as: 


m 


X  +  Y 

where  X,  Y,  and  Z  are  functions  of  mole  fraction,  molecular 
weight  of  the  component,  and  the  collision  integrals  for 
the  unlike-interaction.   Expressions  for  X,  Y,  and  Z  have 
been  given  elsewhere . ' 32> 

In  Table  9  the  calculated  viscosities  and  the  devia- 
tions from  experimental  data  are  presented.   The  average 
deviation  of  calculated  values  from  experiment  is  about 
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1.3  #  which  favorably  compares  with  the  results  of  Mason 

and  Monchick.^30''   The  agreement  is  within  experimental  error 

and  not  worse  than  those  of  the  12-6  potential' 32'    or  modified 

n        i  •    i_         j   -i  (38,  39J 

Buckingham  model        for  nonpolar  spherical  molecules. 

The  diffusion  coefficient  of  binary  mixtures,  Di2, 
is  calculated  by:^32' 

0.002628  (l3/2m)x^z 
D"=  Poa   B.(lfl)  (45) 

12      12 

where  m  is  reduced  molecular  weight,  P  the  pressure  in 
atmosphere,  T  in  °K,  and  a  the  separation  at  zero  potential 

o 

in  A.   The  calculated  values  are  compared  with  experiment  in 
Table  10.   The  average  deviation  of  calculated  values  from 
experimental  data  is  about  5.5$  which  again  confirms  that 
this  model  is  as  good  as  either  the  12-6  potential  or  the 
modified  Buckingham  model  for  nonpolar  molecules. 

It  is  interesting  to  compare  the  parameters  for 
H2O  +  O2  system  with  those  obtained  by  others.   The  values 

o 

by  Eqs.  (31 )  and  (32)  are   012  =  3.^9  A  and   £l2*/k  =  167. 7°K 
which  lie  between  Mason  and  Monchick's^ 3oy  (  °Xz   =  3.072  A, 
e  i2*/k  =  239°K)  for  the  Stockmayer  potential  and  Walker  and 
Westenberg  's  ^4°^  (  o12  =  3.335  A,  E12A  =  8o°K)  for  the 
12-6  potential.   The  systems  involving  H2O  as  one  component 
show  improved  agreement  with  experimental  over  Mason  and 
Monchickfs  results. 
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The  thermal  diffusion  factor  is  a  complex  function  of 
temperature,  concentration,  and  molecular  interaction.   It 
is  hard  to  predict  this  property  theoretically.   Two 
approximations  are  those  of  Chapman-Cowling^33'  and  Kihara- 
Mason.^8'41'   Kihara's  first  approximation  was  used  because 
of  its  simplicity.   This  Is  known  to  give  more  accurate 
results  than  Chapman- Cowling's  first  approximation  and 
differs  from  the  second  approximation  by  a  few  percent. 
In  Kihara's  first  approximation,  the  thermal  diffusion 
factor,  a  _,  is  written  as: 

«   =  (6Ci2*-5)  XlSl  "  X2S2 (46) 

Xl  Ql  +  X2  Q2  +  XlX2  Ql2 

where  x  is  mole  fraction  of  the  component.   Expressions 
for  Si,  S2,  Qi,  Q2  and  Q12  have  been  given  by  Mason^42' 
but  they  are  omitted  here  because  of  their  length.   C12* 
is  the  ratio  of  collision  integrals,  Q*^1'2'    over  ft*^1'1  , 
for  the  unlike-interaction. 

The  experimental  data  for  the  systems  H2O  +  H2,4 
S02  +  H2,   4'  Ar  +  HCl/45'  and  Ne  +  NH3^46^  are  compared 
with  calculated  values  in  Figures  4,  5,  and  6.   Curves  refer 
to  the  calculated  values  and  circles  refer  to  the  experimental 
values.   In  Figure  4  open  circles  and  curve  1  refer  to  values 
at  566°K  and  solid  circles  and  curve  2  to  those  at  368°K. 
In  Figure  6  curve  1  and  open  circles  refer  to  the  values 
for  Ar  +  HC1.   The  temperature  for  S02  +  H2  is  376°K  and  that 
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for  Ar  +  HC1  and  Ne  +  NH3  is  383°K.   The  overall  agreement 
is  fair  and  compares  with  other  calculations  for  polar  gases. 
For  the  H2O  +  H2  system  the  data  of  Whalley  seem  too  uncertain 
with  the  trend  of  the  calculated  values.   However,  the 
Ne  +  NH3  system  shows  a  good  agreement  and  the  inversion 
point  is  correctly  predicted  as  shown  in  Figure  6. 

It  can  be  said  from  the  results  that  the  12-6-6 
potential  can  be  used  for  polar  gas  mixtures  with  success. 
The  12-6-6  model  gives  correlation  at  least  as  good  as, 
or  better  than  the  more  elaborate  model  of  Monchick  and 
Mason  for  mixtures  containing  polar  gases.   When  comparing 
the  result  on  viscosities  of  binary  mixtures  between  these 
two  models,  one  should  bear  in  mind  that  Monchick  and  Mason 
used  the  experimental  data  for  pure  components  to  force 
the  calculated  values  to  fit  at  both  ends  of  the  composition 
range.   When  the  calculated  viscosity  was  used  for  one 
component  in  H2  +  (02^)20   system  the  deviation  was  the 
largest  among  all  the  systems  used  in  their  calculation. 

In  the  case  of  those  systems  containing  NH3  as  one 
component,  the  deviation  Is  the  largest  at  the  polar  end 
of  composition.  This  arises  from  the  fact  that  when  the 
parameters  for  NH3  were  obtained  from  the  viscosities  at 
two  different  temperatures,  the  deviations  at  other  tempera- 
tures were  rather  large.  Thus,  at  low  temperatures  near 
those  used  in  the  fit,  this  trend  disappears. 


As  the  agreement  is  not  very  good  for  the  thermal 
diffusion  factor  which  is  the  most  sensitive  to  the 
unlike-interaction,  it  is  suggested  that  the  combining  rules 
can  be  refined  to  give  better  agreement  between  experiment 
and  calculation.   However,  the  agreement  in  the  diffusion 
and  viscosity  does  not  warrant  any  further  complications 
in  the  combining  rules. 


CHAPTER  V 


SECOND  VIRIAL  AND  JOULE- THOMSON 
COEFFICIENTS  OF  POLAR  GASES 


In  Chapter  III  the  parameters  for  the  angle-averaged 
potential  model  were  determined  from  the  experimental 
viscosities  of  polar  gases.   The  calculations  of  transport 
properties  with  these  parameters  were  very  successful  but 
those  of  second  virial  coefficients  were  not  impressive. 
In  this  chapter  the  parameters  are  determined  from  the 
existing  second  virial  coefficients  of  polar  gases.   The 
validity  of  the  parameters  are  further  tested  by  computing 
the  Joule-Thomson  coefficients  at  zero  pressure  for  several 
molecules  for  which  experimental  data  are  available. 

To  determine  the  parameters,  e  °  and  r°,  from  the 
second  virial  coefficients,  a  similar  method  to  that  of 
Chapter  III  is  chosen.   For  the  12-6  potential  the  second 
virial  coefficient  is  given  by: 


2  2 
B(T)  =  ^ —   *  N  r*3  B*  (47) 


where  the  symbols  have  the  same  meaning  as  before  and  values 
of  B*  have  been  tabulated  elsewhere. ' 33'      Wo  choose  two 
temperatures,  Tx  and  T2,  and  the  corresponding  values  of 
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B(T),  Bi  and  B2,  where  subscripts  1  and  2  refer  to  the 
temperatures.   Eliminating  r*'s  from  these  equations  by 
the  use  of  Eq.  (34)  the  working  equations  of  Eqs.  (35) 
and  (36)  are  obtained  with: 

Ai  =  1.3916xl0+7  V   4(Bi*)2  ( g  x  _  6  2)         (48) 

B2 

l 

and  similarly  for  A2.   The  symbols  are  in  the  same  unit 
as  before.   With  two  known  Bi  and  B2  at  Ti  and  T2,  the 
values  of  e  x*  and   e2*  and  consequently  those  of  rx* 
and  r2*  can  be  evaluated.   Then  by  Eqs.  (38)  and  (39) 
£  °  and  r°  can  be  computed. 

The  potential  parameters  for  several  polar  gases 
determined  by  this  method  are  given  in  Table  11.   Values 
of  dipole  moment  and  polarizability  in  Chapter  III  are  used 
except  the  following. 


m   ,debye 

a  ,A3 

CH3CN 

3.44 

4.47 

CH3CHO 

2.68 

4.43 

C2H5C1 

2.05 

6.4 

The  value  of  a  for  CH3CN  is  taken  from  Moelwyn-Hughes^47' 
and  a  for  CH3CHO  is  calculated  from  the  bond  polarizabilities 
All  other  values  are  given  in  Landolt-Bornstein  Tables. 
The  experimental  second  virial  coefficients  of  S02  are  taken 
from  Kang,  et  al/36^  and  others  from  Rov/linson/ 2G'      The 
values  of  e  1*  and  e  2*  are  found  graphically.   However, 
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they  could  not  be  located  for  CH3OH  and  CH3CI  for  which 
the  curves  by  Eqs.  (35)  and  (36)  ran  parallel. 

The  Joule-Thomson  coefficient  is  defined  as  the 
rate  of  change  of  temperature  with  pressure  in  an  isenthalpic 
expansion.   This  property  is  related  to  the  equation  of 
state  by: 


9  V 
3  T  j 


-   V 


/c 


(49) 


where  "p  is  the  Joule-Thomson  coefficient,  V  is  the  volume, 
and  C  is  the  heat  capacity  of  the  gas.   This  equation  can 
be  expanded  by  the  use  of  the  virial  equation  of  state  and 
in  the  limiting  case  of  p  =  o,  the  Joule-Thomson  coefficient 
at  zero  pressure,  y  °,  is  obtained  as  follows: 


—  o 
y 


CP° 


(T 


dB 
dT 


-  B 


(50) 


where  C  °  is  the  molar  heat  capacity  at  zero  pressure  and 
at  that  particular  temperature  under  consideration.  Since 
the  potential  parameters  are  temperature-dependent  we  have: 


dB 
dT 


dB* 
dT* 


b_(T*  -  r  -B*)  +  TB*  -.  b0 


dbc 

dT 


T*2  dB*  de  *   (51' 
k   dT*  dT 


where  b0  is  (2/3)  nN  o3  and  other  reduced  variables  are 
defined  earlier.   Eqs.  (3l),  (32),  and  (51)  combined  with 
(50)  give  the  Joule-Thomson  coefficient  at  zero  pressure 
for  the  present  model. 
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bo(T*  f|J  -B*)  +  ^^ L«g\« 


6T*     5   dT#| 


(52) 


When  molecules  are  nonpolar,  only  the  first  term  of 

Eq.  (52)  remains  which  is  the  Joule-Thomson  coefficient  of 

nonpolar  molecules.   By  the  knowledge  of  either  the  second 

virial  coefficient  and  its  derivative  with  respect  to 

temperature  or  the  potential  parameters,  the  Joule-Thomson 

coefficient  at  zero  pressure  can  be  easily  calculated. 

Tables  of  T*  ^~     have  been  computed  for  the  12-6  potential/32 

The  Joule-Thomson  coefficients  at  zero  pressure 
calculated  by  using  the  parameters  determined  from  the 
second  virial  coefficients  are  compared  with  the  experiment 
in  Table  12.   The  experimental  data  are  taken  from  Diaz-Pena/48-' 


The  heat  capacities  at  zero  pressure  for  SO2  and  H2O  are 

froir 
(50) 


taken  from  Hougen,  et  al/49'  and  those  of  NH3  from  Osborne, 


et  al. 

The  experimental  and  calculated  second  virial  coeffi- 
cients are  compared  in  Figures  7,    8  and  9  where  the  curves 
are  calculated  values.   The  data  points  used  in  the  determina- 
tion of  the  parameters  are  marked  by  crosses.   For  most  of 
the  molecules  considered  the  agreement  is  within  experimental 
error.   For  NH3,  SO2  and  H2O  the  agreement  is  extremely  good 
but  for  CH3CI  and  (CH3)2C0  the  curvature  of  calculated  curves 
seems  a  little  larger  than  that  of  experiment. 
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The  values  of  e  °/k  obtained  from  the  second  virial 
coefficients  are  always  larger  than  those  from  viscosities. 
As  can  be  expected  r°  values  are  smaller  than  those  of 
Chapter  III  except  H20  for  which  both  e  °/k  and  r°  are 
larger  than  those  obtained  from  viscosities.   In  Table  11 
the  ratios  of  empirical  to  theoretical  dispersion  interaction 
constant  are  also  given.   The  theoretical  values  are  computed 
by  the  theory  of  Mavroyannis  and  Stephen  of  Eq.  (4l).   It 
is  interesting  to  note  that  except  for  H2O  the  ratios  are 
smaller  than  those  of  Chapter  III.   This  trend  is  contrary 
to  the  expectation  based  on  the  fact  that  the  integrals 
of  the  second  virial  coefficient  and  transport  property 
place  emphasis  on  the  attractive  and  repulsive  section  of 
the  potential,  respectively.   However,  the  parameters  are 
consistent  compared  with  those  for  the  Stockmayer  potential. 
In  the  latter  case  a  for  CH3CI  and  CH2C12  are  smaller  than 
that  for  CH4  which  is  most  unlikely. 

The  mean  deviation  in  the  calculations  of  the  Joule- 
Thomson  coefficient  at  zero  pressure  is  10.8$.   For  H20, 
NH3,  and  S02  the  deviations  are  3.3,  16.4  and  16.8$, 
respectively,  which  seems  quite  satisfactory.   The  contribu- 
tion of  the  second  term  in  Eq.  (52)  ranged  from  about  10  $ 
up  to  about  40$  in  this  calculation.   The  agreement  for 
H20  Is  excellent  but  for  NH3  and  S02  the  deviation  seems 
rather  systematic.   The  second  virial  coefficient  and  heat 
capacity  are  reproduced  within  0.5$  with  the  correlations 
used  for  these  molecules.   In  this  respect  the  deviations 
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in  the  Joule-Thomson  coefficient  are  unreasonable.   It 
may  be  possible  that  the  experimental  data,  either  y  or 
B  or  both,  are  in  error.   At  least  they  are  not  completely 
consistent. 

The  angle-averaged  potential  is  as  good  for  polar 
molecules  as  the  12-6  model  for  nonpolar  molecules.   It 
can  reproduce  both  thermodynamic  and  transport  properties 
of  gases  within  experimental  error.   Unfortunately,  it 
apparently  fails  to  correlate  both  properties  with  a  single 
set  of  parameters. 


CHAPTER  VI 


ESTIMATION  OP  12-6-6  POTENTIAL  ENERGY  FUNCTION 

PARAMETERS  FROM  CRITICAL  TEMPERATURES 

AND  POLARIZABILITIES 


In  the  preceding  chapters  the  angle-averaged  pair- 
potential  energy  function  has  been  investigated  to  correlate 
properties  of  polar  gases.   The  model  with  temperature- 
dependent  parameters  proved  to  be  successful.   The  potential 
parameters  are  determined  from  experimental  viscosities 
and  second  virial  coefficients  and  this  is  usually  the  case 
with  other  potentials.   However,  for  any  potential  model 
to  be  truly  useful  in  predicting  properties  of  materials, 
there  must  be  a  method  of  estimating  the  parameters  with 
confidence  from  some  simple,  easily  accessable  properties. 
In  fact,  for  the  Lennard- Jones  12-6  potential  there  have 
been  several  methods  of  estimating  the  parameters,  i.e., 
from  critical  constant,''32''  from  polarizabilities,  ^  51^and 
from  saturated  liquid  densities  at  low  pressure. ^ 52'   In 
this  chapter  the  angle-averaged  potential  parameters  are 
estimated  from  critical  temperatures  and  polarizabilities. 

The  principle  of  corresponding  states  has  been  used 
to  predict  the  potential  parameters  for  12-6  potential 
model.  ^32>      This  principle  is  based  on  several  assumptions4,5  3' 5' 
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among  which  the  following  are  important  to  the  present 
discussion;  (l)  the  molecules  are  spherically  symmetrical, 
(2)  the  potential  energy  is  only  a  function  of  intermolecular 
separation.   Thus,  with  the  12-6  potential  model  which  is 
•good  for  simple,  spherically  symmetrical  molecules,  large 
molecules  show  deviations  from  the  corresponding  states 
theorem  and  should  be  treated  differently/55'56^   As  for 
the  polar  molecules,  the  above  assumptions  cannot  be  true. 
The  present  model,  Eq.  (29),  was  rather  successful  with 
polyatomic,  nonspherical  molecules  as  was  shown  in  earli 
chapters.   Moreover,  the  potential  energy  is  only  a  function 
of  intermolecular  separation  in  the  form  of  Eq .  (30)  although 
it  is  also  an  implicit  function  of  temperature,  dipole 
moments,  and  polarizabilities.   Eq.  (30)  also  satisfies  the 
requirement  of  universal  functional  form  for  the  potential 
energy.   Therefore,  we  can  follow  the  derivation^53' 54)  0f 
the  principle  of  corresponding  states.   Indeed  the  dependence 
of  the  parameters  on  those  properties  does  not  hinder  the 
integration  over  momenta  in  the  partition  function.   From 
the  conditions  for  the  critical  states  it  follows  that  at 
the  critical  state,  kT,/  £*c  and  Vc/r*c3  should  be  universal 
constants.   Here  we  used  the  subscript  c  to  denote  critical 
values. 

For  the  present  purpose  of  estimating  e°  and  r°  only 
e  *c  need  be  known.   Using  the  parameters  in  Chapter  III, 
values  of  E*c  /kTc  are  calculated  and  given  in  Table  13 
along  with  critical  temperatures.   In  view  of  the  nature  of 
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water  it  was  taken  out  in  the  average.   It  was  thus  obtained 
the  relationship, 

e*c/k  =  0.727  Tc  (53) 

The  calculated  values  of  e  *  /k  by  Eq.  (53)  are  compared 
with  experiments  also  in  Table  13.   The  average  deviation 
is  about  16$  .   Since  all  the  molecules  are  large  and 
nonsymmetrical  the  deviation  seems  not  so  bad. 

The  theory  of  Mavroyannis  and  Stephen  for  the  dispersion 
attraction  is  given  in  Eq.  (4l).   This  theory  is  in  the  same 
functional  form  as  that  of  Slater-Kirkwood^ 57 '    but  it  uses 
the  total  number  of  electrons  in  the  atom  rather  than  the 
uncertain  number  of  valence  electron  in  the  latter  theory. 
It  was  shown  in  Chapter  III  that  the  deviation  of  the 
calculated  values  of  e°r°   from  the  empirical  values 


becomes  large  for  large  molecules.   Eq.  (4l)  suggests  that 
when   e°r°   is  plotted  against  a 3Z  on  a  full  logarithmic 
scale,  it  should  yield  a  straight  line  of  slope 0.5.   In 
Figure  10  this  kind  of  plot  was  made.   For  the  polar  molecules 
all  the  information  is  taken  from  the  preceding  chapters. 
The  broken  line,  calculated  by  Eq.  (4l),  does  not  satisfy 
those  empirical  points.   By  the  method  of  least  squares 
an  empirical  equation  is  obtained  for  the  molecules 
considered  here. 

e°r°6  =  h. 9S8xlO"12  (  a3z)°*5985  erg  A6   (5*0 
The  empirical  line  is  also  shown  in  Figure  10. 
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The  parameters,  e  °  and  r°,  can  now  be  estimated. 


■?, 


From  Eqs.  (53)  and  (54)  r*cA  and   e°r°"  can  be  calculated. 
Then  by  Eq.  (3l)   e°/k  can  be  evaluated  with  known  dipole 
moment  and  polarizability.   The  value  of  r°  follows  directly 
from  the  value  of   e°r°6.   In  this  way  the  values  of  e ° 
and  r°  are  calculated  and  compared  with  empirical  values 
in  Table  14.   The  average  deviations  of  the  calculated 
values  of  e °/k  and  r°  are  14.0$  and  5.7$,  respectively, 
which  seem  quite  satisfactory  and  compare  favorably  with 
those  of  16.9$  and  ~$A%  of  Brandt.  ^51^ 

To  check  the  applicability  of  Eqs.  (53)  and  (54)  the 
second  virial  coefficients^58'  and  viscosities     of 
Freons  are  calculated.   The  estimated  parameters  are  given 
in  the  lower  part  of  Table  14.   Dipole  moments  and  polariza- 
bilities  are  taken  from  Landolt-Bornstein  Tables.  3     The 
values  of  dipole  moments  in  debye  unit  are  0.45,  0.51  and 

o 

1.29  and  polarizabilities  in  A3  of  8.24,  6.34,  and  4.50 
for  Freon-11,  Freon-12,  and  Freon-21,  respectively.   The 
polarizability  of  Freon-21,  was  calculated  by  Brandt^51' 
from  the  bond  polarizabilities. 

The  calculated  second  virial  coefficients  and  viscosities 
with  the  estimated  parameters  are  compared  with  experiments 
in  Tables  15  and  16.   The  agreement  is  surprisingly  good 
for  both  properties.   However,  the  overall  agreement  In 
Freon-21  calculations  is  only  fair.   Since  the  polarizability 
comes  in  every  step  of  the  calculation,  we  would  rather 
suspect  the  accuracy  of  the  calculated  polarizability  of  Freon-21. 
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As  for  Eq.  (53)  the  coefficient  of  0.727  is  smaller 
than  that  of  nonpolar  molecules.   According  to  Cook  and 
Rowlinson' s^ 56'  treatments  of  elliptical  and  polar  molecules, 
the  effect  of   asymmetry  and  dipole  moment  makes  the  critical 
temperature  higher  than  that  corresponding  to  spherical, 
nonpolar  molecules.   For  polar  molecules  they  obtained 
the  relationship: 

Td  =  T  °(1  +  4  6  )  (55) 

c     c  c  x 

where  subscript  c  denotes  critical  values  and  superscripts 
d  and  o  denote  dipolar  and  hypothetical  nonpolar  molecule  . 

They  obtained  values  of  6   from  the  consideration  of  the 

c 

corresponding  states  principle  on  vapor  pressure  as 

follows:   0.012  for  H2S,  0.028  for  HC1  and  0.012  for  HBr. 

When  averaged  we  have   6   =  0.0173  which  in  turn  gives 

T  =  1.07  T  °.   From  Eq.  (53)  and  that  for  nonpolar  gases, 

i.e.,   e*  /k  =  O.78  T  as  used  by  them,  we  obtain 

T  =  1.07  T  °  which  in  some  sense  confirms  the  validity 
c         c 

of  Eq.  (53). 

6 

Turning  to  the  empirical  equation  for  e°r°  ,  there 
can  be  some  explanation*51'  on  the  exponent.   That  is, 
the  omission  of  higher  dipole-dipole  interaction  terms 
should  be  counterbalanced  by  the  use  of  higher  exponent. 
The  value  of  0.5985  lies  between  theoretical  value  of  0.5 
and  Brandt's  value  of  0.62.   However,  Brandt  used  Slater- 
Kirkwood's  formula  instead  of  Mavroyannis-Stephen.   In 
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the  present  potential  model  dipole-induced  dipole  interaction 
is  accounted  for  and  the  exponent  seems  reasonable. 

In  conclusion  it  is  felt  that  the  method  of 
estimating  parameters  from  critical  temperatures  and 
polarizabilities  seems  satisfactory.   Nevertheless,  if 
there  are  empirical  parameters  available,  they  should  be 
preferred  to  those  estimated. 


CHAPTER  VII 


CORRESPONDING  STATES  CORRELATION  FOR 
TRANSPORT  PROPERTIES  OP  GASES 


The  principle  of  corresponding  states  has  been 
widely  used  for  the  correlation  of  both  equilibrium  and 

Up)  ( ■  ■    ) 

nonequllibrium  properties  of  fluids.      Pitzerv  3'  has 
proved  the  principle  on  the  statistical  mechanical  basj 
by  moulding  the  partition  function  Into  a  reduced  variable 
form  with  several  assumptions.   The  reducing  parameters 
can  be  the  critical  constants,  the  potential  parameters 
or  something  else.   When  the  Lennard- Jones  12-6  potential 
is  chosen  as  the  universal  potential  energy  function, 
the  simple  spherical  nonpolar  molecules  follow  exactly 
the  principle  of  corresponding  states.   For  polyatomic 
molecules  Trappeniers     formulated  a  new  principle  of 
corresponding  states  by  Improving  the  crude  spherical 
approximation  for  the  potential  field.   Cook  and  Rowlinson^ " 
treated  the  elliptical  and  polar  molecules  by  the  use  of 
angle-dependent  potential  model.   Guggenheim  and  Wormald 
accounted  for  the  systematic  deviations  of  hydrocarbons 
from  the  principle  of  corresponding  states  by  introducing 
an  additional  parameter.   Thodos  and  co-workers  *"  G1'  e2' 
derived  equations  for  transport  properties  of  gases  by 
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dimensional  analysis  with  the  critical  constants.   Their 
results  were  successful  for  simple  molecules. 

According  to  the  principle  of  corresponding  states 
the  reduced  properties  of  materials  should  be  a  universal 
function  of  reduced  temperature  and  pressure  for  the 
classical  mechanics.   With  the  potential  energy  function 
parameters  as  reducing  variables  we  have  the  reduced 
viscosity  and  diffusion  coefficients  as  follows, ' 63' 

a*  =   n  03/M  "2"  e  2  (56) 

D*  =  DM  k.   /  a    e  2  (57) 

where  the  same  symbolism  as  before  is  used.   The  main 
drawback  in  the  above  scheme  is  that  the  potential 
parameters  should  be  known  for  the  molecules  under  considera- 
tion.  Unfortunately,  potential  parameters  are  not  always 
known.   Thus,  instead  of  them  the  critical  constants  are 
frequently  used. 

In  this  chapter  the  saturated  liquid  density  parameters 
at  low  pressure  are  used  as  the  reducing  variables.  Reed^64' 
has  shown  that  the  dependence  of  saturated  liquid  density 


on  temperature  for  T  <0.8t_  can  be  expressed  as: 


(58) 


where  pis  the  molar  density,  A  and  B  are  positive  consta 

(   R?] 

Reed  and  McKinleyv  =<=:/  have  shown  that  these  parameters  are 
related  to  the  parameters  of  the  generalized  n-m  potential 
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by  employing  the  smoothed  potential  cell  model  for  liquid 
phase  and  n-m  potential  for  gas  phase.   Thus,  the  relation- 
ships are: 


A 

—  =  (  eA)  Ci  (59) 

B 


A  =  C2/r*6  (60) 

where  C's  are  constants  characteristic  to  the  exponents 
n  and  m.  By  combining  Eqs.  (56)  and  (57)  with  (59)  and 
(60)  the  reduced  properties  of  gases  are  obtained. 


^         1  ,,      i    5//<3   -i 
n    =  n  1  MM  2  A     k  ^)       (61) 

^  ^       i   3/2  ,,   11/6   3/2 

D  ^  =  DM  2  B   /(A      k     )       (62) 

The  reduced  viscosity  and  diffusion  coefficients  are  now 
functions  of  reduced  temperature  and  pressure.   The 
reduced  temperature,  T  ^  is  here  TB/A. 

The  experimental  viscosities  of  simple  molecules 
are  reduced  by  Eq.  (61).   These  are  at  pressure  one 
atmosphere  or  less.   The  values  of  A  and  B  used  are  given 
in  Table  17  where  the  sources  of  experimental  density  data 
are  also  given.   By  the  least-squares  method  the  reduced 
viscosities  are  obtained  in  terms  of  reduced  temperatures 
for  eight  molecules. 
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n   =  91.6784  +  1145. 764t  ^  -  91.2258T  ^ 

+  4.1359T  ^  3  -  0.00323T  ^  5  (63) 

This  equation  is  obtained  for  the  temperature  range  of 
0.4^  T  ^  ^18.4.   Noting  that  reduced  viscosities  of 

Ar,  Kr,   02  and  N2  fall  close  to  each  other  a  different 
equation  is  obtained  for  these  molecules, 

V  =  -9.0278  +  1348.23877T  ^  -  183.03746T  ^  2 

+  14.47759T  ^3  -  0.04108T  ^  5  (64) 

for  the  temperature  range  of.0.4^T  r  ^  9.1. 

The  experimental  diffusion  coefficients  at  one 
atmosphere  are  also  reduced  by  Eq.  (62).   When  the 
reduced  diffusion  coefficients  are  plotted  against 
reduced  temperatures  it  appears  that  they  could  be  best 
represented  by  two  straight  lines.   By  the  least-squares 
fit  the  following  two  equations  are  obtained. 

D  ^  x  IO20  =  9.12T  *   1,9°39  ,  0.3  s  ^   i  5.5     (65) 
D  ^  x  102°  =  11.56T  ^  *-e*2^  3.5  ±  T  ^  49.3    (66) 

The  results  of  viscosity  calculations  are  given  in 
Table  18  where  the  number  of  data  points,  temperature  range, 
percent  deviation  and  the  sources  of  experimental  data 
are  also  given.   For  the  168  data  points  of  eight  molecules 
used  in  the  calculation  the  average  deviation  is  3.65   . 
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In  general  the  deviation  becomes  large  as  the  molecules 
become  more  complex.   The  reduced  viscosities       and 
Ar  at  low  temperatures  are  very  close  but  at  nig]      ma- 
tures the  Ar  viscosities  fall  below  those  of  Ne .   This 
trend  has  been  noticed  previously  by  others.  ^  GG}  r'T '      In 
this  calculation  the  reduced  viscosities  of  N2  and  O2  are 
very  close  to  those  of  Ar.   Therefore,  the  experimental 
data  of  Ne  at  high  temperatures  are  rather  suspected  in 
spite  of  the  conclusion  of  Mason  and  RIceVOD/  that  the 
argon  viscosities  are  low  at  high  temperatures  in  view 
of  the  results  of  energy  beam  scattering  experiments.   The 
discrepancies  in  the  reduced  viscosities  of  Kr  and  Xe 
were  also  noted  by  them. 

When  only  Ar,  Kr,  02  and  N2  are  considered  the 
average  deviation  is  l.l8$  for  105  data  points.   Here 
Kr  shows  the  largest  deviation  among  them  all. 

The  summary  of  diffusion  coefficient  calculations 
is  given  in  Table  19.   For  51  data  points  of  seven  molecules 
the  overall  deviation  is  K.2%   .   The  agreement  is  good 
for  simple  molecules  such  as  Ar,  Kr,  02  and  N2 .   However, 
it  is  poor  for  CH4  and  CF4.   It  should  be  noted  that  the 
scatter  of  the  experimental  data  is  quite  large  for  most 
of  the  molecules  considered  here.   Thus,  the  deviation  seems 
rather  satisfactory. 
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In  this  correlation  it  has  been  implicit  that  all  the 
molecules  obey  a  particular  n-m  potential  energy  function. 
Yet  it  is  a  well-known  fact  that  this  is  not  true. 
That  is,  the  12-6  potential  is  good  for  small  spherical 
molecules  while  28-7  potential  is  better  for  polyatr 
molecules  than  the  12-6  model.   Thus,  from  the  theoretical 
consideration  it  is  hard  to  expect  that  the  molecules 
considered  here  should  follow  the  same  corresponding  st 
As  illustrated  at  the  beginning  of  this  chapter,:- 
molecules  should  be  treated  separately  and  we  should 
be  content  with  the  results  of  these  calculations  which 
compare  favorably  with  others. ^ ex' 61' 
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Table  2 
Potential  Parameters  for  Nonpolar  Molecules 


Molecule 

Potential 

E  /k,°K 

a,  A 

c 

Re: 

Ar 

L-J 

117.7 

3.504 

- 

16 

Morse 

119.53 

3.5149 

4.533 

Kr 

L-J 

164.0 

3.827 

- 

16 

Morse 

192.24 

3.6257 

4.848 

79 

C02 

L-J 

198.2 

4.328 

- 

'i 

SP6 

L-J 

I88.7 

5.910 

- 

32 

CF4 

L-J 

151.5 

4.744 

- 

16 

CH4 

L-J 

148.9 

3.783 

- 

16 

Morse 

203.5 

3.6463 

5.374 

79 

C2H6 

L-J 

230.O 

4.418 

- 

32 

Morse 

469.01 

3.7417 

6.721 

17 

C3H8 

L-J 

254.0 

5.061 

- 

32 

Morse 

527.81 

4.3212 

6.496 

17 

C4H10 

L-J 

410.0 

4.997 

- 

32 

Morse 

650.0 

5.0665 

7.751 

17 

NeoC5Hi2 

L-J 

232.5 

7.445 

- 

16 
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Table  3 
Critical  Volumes  and  Polarizabilities 


V  , cc./mole 
c 

Ref . 

a  ,A3 

a  Re; 

Ar 

74.52 

80 

1.626 

34 

Kr 

92.29 

80 

2.456 

34 

CO2 

94.2 

25 

2.65 

34 

SF6 

199 

25 

6.55 

34 

CF4 

138 

64 

4.02 

34 

CH4 

99 

25 

2.60 

34 

C2H6 

148.2 

25 

4.47 

34 

C3H8 

200.6 

25 

6.29 

34 

C4H10 

255.1 

25 

8.12 

34 

NeoC5Hia 

303.5 

25 

10.36 

34 
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Table  6 
Second  Virial  Coefficient  for  Polar  Gases 


NH3 


-B  cc/mole 

t°c. 

calc  . 

expt . 

-30 

329 

560 

0 

232 

367 

30 

VJk 

261 

60 

136 

197 

100 

101 

143 

150 

73 

103 

200 

54 

77 

250 

40 

59 

H20 


-B 

cc 

./mole 

t°c. 

calc . 

expt . 

150 

534 

284 

200 

315 

197 

300 

154 

112 

400 

92 

72 

SO; 


CH3CI 
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Table  6  (Continued) 

-B  cc 

,/mc 

ile 

t°c. 

calc  . 

expt . 

10 

328 

500 

20 

304 

452 

30 

283 

kok 

40 

265 

367.5 

50 

248 

332.8 

75 

212 

279 

100 

184 

232.5 

125 

160 

201.0 

150 

141 

171.1 

175 

124 

144.1 

200 

110 

125.8 

-B  cc 

./mole 

t°c. 

calc  . 

expt . 

-34 

508 

764 

-23 

457 

668 

-18 

436 

637 

+10 

344 

500 

38 

279 

401 

65 

232 

320 

93 

194 

265 

121 

164 

214 

149 

140 

184 

177 

120 

155 
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Table  7 


Empirical  and  Theoretical  Dispersion 
Attraction  Constants 


erg.  x  A6  x  10: 


Polar 
Molecules 

Empirical 

Equation  12 

Empirical 
Theoretical 

NH3 

74.7 

67.6 

1.105 

H20 

45.7 

36.2 

1.26 

S02 

383 

255 

1.50 

HC1 

127 

114 

1.11 

HI 

427 

587 

0.727 

CH3CI 

468 

312 

1.50 

CH2CI2 

1295 

672 

1.93 

CHCI3 

2160 

1131 

1.91 

CH3OH 

294 

155 

I.89 

C2H5OH 

843 

427 

1.97 

(CH3)2C0 

1436 

566 

2.54 

Nonpolar 
Molecules 

Ne 

4.50 

4.38 

1.03 

Ar 

54.6 

55.36 

0.99 

Kr 

116 

145.3 

0.80 

Xe 

281 

369.9 

0.76 

N2 

62.8 

53.75 

1.17 

o2 

51.1 

49.09 

1.04 

CH4 

128.9 

88.2 

1.46 

CCI4 

3731 

2004 

1.86 
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Table  8 

Potential  Parameters  and  Molecular  Constants 
for  Nonpolar  Molecules 


Molecules     e/k.,°K  r,A  a  ,A3  Ref 

H2  57-0  5.287  0.79  50 

N2  91.5  4.132  1.76  50 

02  115.0-  5.855  1.60  50 

GH4  1^4,0  4,261  2.60  50 

C2H4  205.0  4.750  4.26  50 

C02         •        190.0  4.485  2.65  50 

CCI4  527.0  6.601  10.50  50 

He  10.22  2.869  0.204  50 

Ne  55.7  5.151  0.595  50 

Ar  124.0  5.857  1.626  50 
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Table  9 
Viscosities  of  Binary  Mixtures 


System 

T°K 

X  polar* 

nxlO7 

calc . -expt . 
%   dev. 

Ref . 

NH3+H2 

293. 16 

0.9005 

1003 

-0.1 

86 

O.7087 

1045 

-0.2 

0.5177 

1080 

0.0 

0.2975 

1088 

+0,1 

0.2239 

1073 

+0.1 

0.1082 

1013 

+0.2 

373.16 

0.9005 

1335 

+2.8 

86 

0.7087 

1361 

+2.1 

0.5177 

1372 

+1.3 

0.2975 

1338 

+0.7 

0.2239 

1303 

+0.3 

0.1082 

1207 

+0.2  - 

473.16 

0.9005 

1726 

+4.0 

86 

.  0.7087 

1730 

+3.0 

0.5177 

1711 

+2.1 

0.2975 

1627 

+1.1 

0.2239 

1572 

+0.8 

0.1082 

1438 

+0.4 

523. 16 

0.9005 

1906 

+4.4 

86 

0.7087 

1900 

+3.43 

0.5177 

1865 

+2.3 

0.2975 

1756 

+1.1 

0.2239 

1690 

+0.7 

X  =  mole  fraction 
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Table  9  (Continued) 


System 

T°K 

X  polar 

n  xlOT 

calc . -expt . 
$dev. 

Ref 

NH3+N2 

295. 16 

O.8885 

1094 

+0.2 

86 

0.7147 

1261 

+0.6 

0.5658 

1392 

+0.7  ■ 

0.2920 

1594 

+0.6 

0.1111 

I698 

+0.5 

575.16 

0.8885 

1434 

+2.6 

86 

0.7147 

1601 

+2.0 

0.5658 

1731 

+1.2 

O.2920 

1950 

+0.5 

0.1111 

2055 

+0.1 

475.16 

0.8885 

1830 

+3.5 

0.7147 

1990 

+2.3 

0.5658 

2115 

+1.4 

0.2920 

230  6 

+0.4 

0.1111 

2409 

0.0 

525.16 

O.8883 

2013 

+3.8 

86 

0.7147 

2169 

+2.7 

0.5638 

2291 

+1.8 

0.2920 

2480 

+0.8 

0.1111 

2583 

+0.4 

NH3+02 

295.16 

0.8755 

1159 

-0.3 

86 

0.7079 

1545 

-0.5 

- 

0.W6 

1602 

-0.1 

System  T°K  X  polar  nxl07      $dev.  Ref 
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Table 

9  (Continued) 

X  polar 

nxl07 

calc . -expt . 
#dev. 

O.2986 

1783 

0.0 

0.1351 

1927 

+0.2 

16 

0.8755 

1489 

+2.1 

0.7079 

1707 

+1.0 

0.4786 

1982 

+0.5 

O.2986 

2174 

+0.2 

0.1351 

2329 

+0.1 

16 

0.8755 

1896 

+3.0 

0.7079 

2123 

+1.8 

0.4786 

2409 

+0.8 

O.2986 

2611 

+0.3 

0.1351 

2778 

+0.2 

66   • 

0.9 

990 

-1.8 

0.8 

1019 

-1.9 

0.7 

1043 

-1.7 

0.6 

1063 

-1.3 

0.5 

1078 

-1.2 

0.4 

1087 

-1.1 

0.3 

1092 

-1.2 

0.2 

1090 

-1.3 

0.1 

1083 

-1.5 

86 


86 


NH3+CH4  287.66       -0.9  990  -1.8  100 
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Table  9  (Continued; 


System 

T°K 

X  polar 

n  xlO7 

calc . -expt . 
%  dev. 

Ref 

NH3+C2H4 

295.16 

0.8867 

1013 

+1.2 

86 

0.8071 

1030 

+1.7 

0.6961 

1047 

+2.4 

0.5172 

10  5^ 

+2.7 

0.2993 

1049 

+2.1 

0.1096 

1027 

+1.2 

373.16 

0.8867 

1335 

+3.2 

86 

0.8071 

1342 

+3.2 

0.6961 

1345 

+3.1 

0.5172 

1338 

+2.7 

0.2993 

1312 

+1.6 

0.1096 

1277 

+0.6 

86 

473.16 

0.8867 

1712 

+3.9 

O.807I 

1704 

+3.4 

0.6961 

I689 

+3.1 

0.5172 

1657 

+2.2 

0.2993 

1609 

+0.9 

0.1096 

1561 

+0.0 

523.16 

0.8867 

1886 

+4.3 

86 

0.8071 

1871 

t-3.7 

0.6961 

1848 

+3.2 

0.5172 

1805 

+2.3 

0.2993 

1747 

+1.0 

0.1096 

1693 

+0.2 
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Table   9    (Continued) 


System 

T°K 

X  polar 

n  xlO7 

calc . -expt . 
%    dev. 

Ref 

HCI+H2 

29^.6 

0.8220 

1451 

-0.7 

94 

0.7179 

1459 

-0.6 

0.5042 

1460 

-0.7 

0.2031 

1331 

-0.8 

327. 16 

0.8220 

1617 

-0.6 

94 

0.7179 

1623 

-0.6 

0.5042 

1614 

-0.7 

0.2031 

1451 

-1.4 

372.16 

0.8220 

1839 

-0.5 

9^ 

0.7179 

1840 

-0.8 

0.5042 

1816 

-0.8 

0.2031 

1607 

-1.4 

94 

427-16 

0.8417 

2099 

0.0 

O.6989 

2094 

-0.5 

0.5092 

2057 

+0.2 

0.2409 

1855 

-0.6 

473.16 

0.8417 

2307 

-0.2 

9^ 

0 . 6989 

2296 

-0.3 

0.5092 

2245 

-0.7 

0.2409 

2008 

-0.8 

523. 16 

0.79^7 

2520 

-0.3 

9^ 

0.6312 

2489 

-0.7 

0.5178 

2446 

-0.3 

0.2991 

2256 

-1.1 

Table 

9  (Continued) 

System 

T°X 

X  polar 

n.  xlOy 

calc . -expt . 
%    dev. 

Ref 

HCI+CO2 

291. 16 

0.9 

1426 

-2.3 

100 

0.8 

1437 

-2.4 

0.7 

1445 

-2.6 

0.6 

1451 

-2.7 

0.5 

1456 

-2.9 

0.4 

1459 

-2.9 

0.3 

1460 

-2.9 

0.2 

1459 

-2.7 

0.1 

1457 

-2.5 

SO2+H2 

290. 16 

0.8215 

1273 

-1.5 

90 

0.5075 

1333 

-1.3 

0.2963 

1346 

-1.7 

0.2286 

1330 

-1.0 

0.1676 

1294 

-0.8 

318. 16 

0.8028 

1409 

-1.1 

90 

0.5075 

1462 

-0.9 

0.2963 

1464 

-1.3 

0.2286 

1441 

-0.8 

0.1676 

1397 

-0.9 

343.16 

0.8028 

1527 

-0.5 

90 

O.6999 

1545 

-0.8 

0.6175 

1559 

-1.0 

0.4823 

1577 

-0.6 
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Table  9  (Continued) 

System      T°K       X  polar    nxlO7  %  dev.         Ref 


365.16    0.8028    1628         -0.3      90 


X  polar 

nxlO7 

calc . -expt 
%  dev. 

0.2963 

1568 

-1.8 

0.2306 

1540 

-0.7 

0.1676 

1487 

-0.9 

0.1657 

1485 

-1.3 

0.8028 

1628 

-0.3 

O.6999 

1646 

-0.1 

0.6175 

1659 

-1.0 

0.4823 

1674 

-0.5 

0.2306 

1624 

-1.0 

0.1676 

1564 

-0.6 

0.1657 

1562 

-1.0 

0.6760 

1793 

-0.2 

0.4698 

1813 

-0.1 

0.3265 

1796 

-0.3 

0.1636 

1673 

-0.7 

0.6760 

1946 

+0.2 

0.4698 

I960 

0.0 

0.3265 

1934 

-0.4 

0.1676 

1797 

-0.3 

0.1512 

1769 

+1.2 

0.6760 

2115 

-0.1 

0.4905 

2123 

+0.1 

0.3265 

2085 

-0.6 

0.1512 

1893 

-3.1 

397.16    0.6760    1793         -0.2      90 


432. 16     0.6760     1946         +0.2       90 


472. 16     0.6760     2115         -0.1       90 
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Table  9  (Continued) 

calc . -expt . 
System      T°K       X  Polar    nxlOY  <?o  dev .  Ref, 

SO2+CO2     288.96    0.9       1262         -2.0      100 


CH2Cla+CCl4  293.15   0.8425    987         -3.3      101 


0.9 

1262 

-2.0 

0.8 

1289 

-2.1 

0.7 

1316 

-1.6 

0.6 

1340 

-1.7 

0.5 

1364 

-1.4 

0.4 

1385 

-1.6 

0.3 

1404 

-1.7 

0.2 

1421 

-1.8 

0.1 

1434 

-2.0 

0.8425 

987 

-3.3 

O.6985 

989 

-2.7 

0.5014 

988 

-2.5 

0.3114 

985 

-1.5 

0.1484 

981 

-1.0 

0.6739 

1196 

-1.3 

0.3649 

1190 

-0.1 

0.8485 

1403 

-1.5 

0.7H8 

1401 

-0.1 

0.5264 

1396 

-1.1 

0.2904 

1387 

+0.4 

0.1261 

1380 

+0.9 

353.16   0.6739    1196        -1.3     101 


413.43   0.8485    1403        -1.5      101 


79 

Table  10 
Diffusion  Coefficients  of  Binary  Mixtures 

calc . -expt . 
System         T°K         D12    %  dev.  Ref , 

H2O+H2         293.11  0.828  -2.6         102 

322.66  O.981  -3.0 

365.56  1.221  -1.0 

372.50  1.261  -1.6 

307.26  0.900  -11.8        103 

328.56  1.013  -9.6 

352.66  1.147  -4.4 

307. 16  0.900  -1.6         104 

328.66  1.014  +5.5 

H2O+N2         307.56  0.254  -0.8         103 

328.56  O.287  -5.3 

352.16  0.325  -9.2 

328.96  0.288  -8.0         104 

349. 16  0.321  -9.3 

H2O+O2        308. 06     0.258        -8.5        103 

-8.2 
-6.0 

-13.1        40 
-8.7 
-6.8 
-6.6 
.   -6.7 


308. 06 

0.258 

328.96 

0.292 

352.36 

0.331 

400.00 

0.417 

500.00 

0.621 

700.00 

1.118 

900.00 

1.718 

1050.00 

2.230 
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Table  10 

( Continued) 

calc . -expt 

System 

T°K 

D12    % 

dev. 

H2O+CO2 

296.11 

0.175 

+6.7 

365.56 

0.261 

+5.2 

372.54 

0.270 

+4.2 

307.46 

0.188 

-6.9 

328.56 

0.213 

+0.9 

352.36 

0.243 

-0.8 

328.66 

0.213 

+7.6 

H2O+CH4 

307.66 

0.257 

-12.1 

328.76 

0.291 

-12.1 

352.26 

0.331 

-  7.0 

H2O+C2H4 

307.76 

0.186 

-  8.8 

328.46 

0.211 

-  9.4 

352.56 

0.241 

-  2.4 

H20+He 

307. 16 

0.880 

-  2.4 

328. 16 

0.984 

-  2.7 

352.46 

1.109 

-  1.1 

298. 16 

0.837 

-  7.8 

NH3+H2 

298. 16 

0.785 

+  0.3 

328.16 

0.927 

-  1.7 

358. 16 

I.O78 

-  1.4 

263.0 

0.630 

+10.5 

328.0 

0.926 

+  2.9 

Ref , 
102 

103 

104 
103 


103 


103 

105 
30 

106 


8l 


Table  10  (Continued) 


System 


NH3+N2 


SO2+H2 


SO2+N2 
SO2+CO2 


C2H5OH+H2 
C2H5OH+CO2 
C2H50H+He 
(CH3)2C0+H2 


T°K 

Dia 

calc . -expt 
fodev. 

395.0 

1.276 

-   5.5 

473.0 

1.734 

-  6.3 

298. 16 

0.229 

-  4.3 

328. 16 

0.272 

-  4.6 

358. 16 

0.318 

-  3.0 

285.56 

0.528 

+  0.6 

263.00 

0.456 

+  6.0 

323.00 

0.657 

+  7.7 

473.00 

I.272 

+  3.4 

263.00 

263.00 
343.00 
473.00 

340..  5 
340.16 
298. 16 
296. 16 


0.101 

0.069 
0.116 
0.212 

0.623 
0.107 
0.459 
0.411 


-  2.9 


+ 

7, 

,8 

+ 

7, 

,4 

+ 

8, 

.7 

+ 

6, 

,3 

+ 

1, 

.0 

- 

7, 

,1 

_ 

3 

.1 

Ref , 


30 


107 
106 


106 


106 


102 
102 
105 
102 
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Table  11 

Potential  Parameters  for  the  12-6-6  Potential 
Determined  rom  the  Second  Virial  Coefficients 

e  °r°°  Emp.. 
e  °/k,°K     r°,A      e  °r°S  Theor 

NH3  203.4  3.601  0.906 

H20  234.8  3.670  2.188 

S02  267.1  3.843  0.456 

CH3C1  248.0  4.510  0.923 

C2H5CI  344.7  4.785  0.9S2 

CH3CN  157.7  5.447  2.039 

CH3CHO  167.9  5.040  1.322 

(CK3)2C0  284.1  5.216  1.395 
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Table  12 
Joule-Thomson  Coefficients  at  Zero  Pressure 


System 

T°C. 

Expt . , °K/atm. 

Expt . -Calc . 

Expt  .-Calc 
%   dev. 

NH3 

-18 

4.954 

-1.590 

-32.1 

-13 

4.568 

-1.354 

-29.6 

-  5 

4.045 

-1.054 

-26.0 

10 

3.321 

-0.643 

-19.4 

30 

2.591 

-0.376 

-14.5 

50 

2.091 

-0.221 

-10.6 

70 

1.725 

-0.132 

-  7.6 

110 

1.224 

-0.055 

-  4.5 

150 

0.938 

-0.033 

-  3.5 

S02 

50.7 

3.020 

-0.158 

-  5.2 

101.6 

2.166 

+0.270 

+12.5 

150.5 

1.645 

+0.350 

+21.3 

200.2 

1.315 

+0.372 

+28.3 

H20 

100 

6.37 

-0.52 

-  8.2 

150 

3.85 

+0.15 

+  3.9 

200 

2.48 

+0 .  14 

+  5.6 

250 

i.69 

+0.06 

+  3.5 

300 

1.20 

0 

0.0 

350 

0.91 

-0.02 

-  2.2 

400 

0.73 

-0.01 

-  1.4 

450 

0.595 

-0.013 

-  2.2 

500 

0.490 

-0.016 

-  3.3 

550 

0.415 

-0.012 

-  2.9 
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Table  13 
Empirical  and  Estimated  Values  of  c*  /kTc 


T  ,  °K 

Emp.  e*c  /k,°K 

e*   /kTc 
c 

Calc .-Emp 
$dev. 

NH3 

405.6 

3^7.8 

0.857 

-15 

H20 

647.0 

855.9 

1.323 

-45 

S02 

430.4 

327.7 

0.761 

-*.5 

(CH3)2CO 

508.7 

247.4 

0.486 

+49 

HC1 

324.6 

313.2 

0.965 

-25 

HI 

424.6 

319.8 

0.75^ 

-3.5 

CH3CI 

416.3 

299.6 

0.720 

-0.9 

CH2CI2 

510.0 

377.^ 

0.7^0 

-1.7 

CHCI3 

536.6 

379.5 

0.707 

+2.8 

CH3OH 

513.2 

335.7 

0.654 

+11 

C2H5OH 

516.3 

320.9 

0.621 

+17 

Average 

0.727 

15.9 
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Table  17 
Liquid  Density  Pars 


AxlO4 

BxlO7 

Mole2  x  cc.~2 

Mole2  x  cc. 

"2  x  °K_1 

Ref . 

Ne 

64.9618 

1088.63 

108 

Ar 

22.2310 

113.60 

108 

Kr 

14.&O06 

54.7944 

108 

Xe 

10.0281 

27.2227 

108 

^2 

22.903 

112.84 

65,  76 

N2 

15.540 

93.779 

65,  76 

CH4 

12.2598 

47.2576 

64 

CF4 

6.4629 

21.4311 

64 
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Table  19 
Reduced  Diffusivities  of  Simple  Gases 

%   Dev.  Ref . 

Ne  1.3-5.9        5.25  121 

Ar         7  0.4-1.8        3.72  121 

Kr         1  l.i  2.59  32 

02        24  0.4-49.3      3.43  121,122 

N2         5  0.5-2.1        3.60  121 

CH4        5  0.3-1.4        4.67  121 

1.0-1.1       9.56  123 

4.20 


No.    of 

Pt£ 

Data 

T  r  Range 

5 

1.3-5.9 

7 

0.4-1.8 

1 

1.1 

24 

0.4-49.3 

5 

0.5-2.1 

5 

0.3-1.4 

4 

1.0-1.1 

51 

0.3-49.3 
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